Draft version February 2, 2008 

Preprint typeset using I^T^ style emulateapj v. 04/03/99 



THE DISCOVERY OF A SPATIALLY-RESOLVED SUPERNOVA REMNANT IN M31 WITH Chandra 
Albert K.H. Kong, Michael R. Garcia, Francis A. Primini and Stephen S. Murray 

Harvard-Smithsonian Center for Astropliysics, 60 Garden Street, Cambridge, MA 02138; alcong@cfa.liarvard.edu 

Draft version February 2, 2008 

ABSTRACT 

Chandra observations of M3 1 allow the first spatially resolved X-ray image of a supernova remnant (SNR) in an 
external spiral galaxy. CXOM31J004327.7H-411829 is a slightly elongated ring-shaped object with a diameter of 

1 1" (42 pc). In addition, the X-ray image hints that the chemical composition of the SNR is spatial dependent. 
The X-ray spectrum of the SNR can be well fitted with a Raymond-Smith model or a non-equilibrium ionization 
model. Depending on the spectral model, the 0.3-7 keV luminosity is between 3.2 x 10^*^ ergs"' and 4.5 x 10^^ 
ergs"'. The age of the SNR is estimated to be 3210-22300 years and the number density of ambient gas is 
~ 0.003-0.3 cm""*. This suggests that the local interstellar medium around the SNR is low. 

Subject headings: galaxies: individual (M31) — supernova remnants — X-rays: ISM 



1. INTRODUCTION 

X-ray emission from supernova remnants (SNRs) in our 
Galaxy has been extensively studied in the past twenty years. 
Pre-Chandra observations allowed many SNRs to be resolved, 
revealing several different morphologies. The detailed study of 
SNR morphology has revealed the interaction of ejected mate- 
rials with the interstellar medium (ISM) in detail, and also has 
provided insights into the evolution of the SNR progenitor stars. 
Chandra observations of Galactic SNRs provide exquisite de- 
tail, allowing the chemical stratification of the progenitor to be 
revealed (e.g., Canizares et al. 2001) and identifying the central 
neutron star and its associated synchrotron nebula (e.g., Slane 
et al. 2000). Even with these new insights, studies of Galactic 
SNRs can be limited by the lack of reliable distance estimates 
and generally high interstellar absorption. These limitations are 
overcome by studies of SNRs in nearby extragalactic systems, 
the closest of which (the Large and Small Magellanic Clouds) 
are near enough to still allow Chandra to reveal the structure in 
detail (see e.g., Hughes 2001). 

Because M31 is the closest {r-^ 800 kpc) galaxy which shares 
similar morphology and size with the Milky Way, it is the best 
candidate for a comparative study of SNR. Optical surveys for 
SNR in M3 1 typically identify candidates based on [S II] and 
Hq, imaging observations, e.g. d'Odorico, Dopita, & Benvenuti 
(1980), Blair, Kirshner, & ChevaUer (1981, 1982), Braun & 
Walterbos (1993) and Magnier et al. (1995). About 200 SNR 
candidates have been identified by these surveys, 27 of which 
have been confirmed via spectroscopic observations (see Blair 
et al. 1981, 1982). In a recent extensive ROSATPSPC survey 
(Supper et al. 2001) of M31, 16 SNRs were identified within 
a total area of ~ 10.7 deg^ by cross-correlating with the these 
optical catalogs. Their X-ray luminosities (0. 1-2.4 keV) range 
from - lO^*" to - 10" ergs"'. 

More recendy, Kong et al. (2002) detected two SNRs 
(CXOM3 1 J004327 .7H-4 1 1 829 and CXOM3 1 J004253 .5h-4 1 2553) 
in the central ~ 17' x 17' region of M31 with Chandra, both 
of which were previously identified by ROSAT. It is worth not- 
ing that CXOM31 J004327.7H-41 1829 was also identified as a 
SNR in the Einstein data (Blair et al. 1981). In a recent Chan- 



dra observation of M3 1 core region, we found that one of the 
SNRs is extended. In this Letter, we report on observations 
of the first SNR in M3 1 to be resolved at X-ray wavelengths, 
CX0M31 J004327.7H-411829. 

2. OBSERVATIONS AND DATA REDUCTION 

CXOM31 J004327.7H-411829 was observed several times 
with Chandra Advanced CCD Imaging Spectrometer (ACIS- 
I) under a M31 monitoring program during 1999-2001 (see 
Kong et al. 2002). CXOM31 J004327.7H-41 1829 is about 8' 
from the aim point of these observations, and therefore the in- 
strumental point spread function (PSF) is substantially poorer 
than if it was on-axis. In a recent 5ks ACIS-I observation of a 
ti-ansient in M31 taken on 2001 August 31 (Kong et al. 2001), 
CX0M31 J004327.7H-41 1829is < 4' from the aim-point, where 
the PSF is sufficiently small to allow us to resolve this ob- 
ject. Additional observations were obtained with the ACIS-S 
on 2001 October 5 with total integration time of 37.7ks; al- 
though CXOM31J004327.7H-411829 is in the ACIS-S4 front- 
illuminated (FI) chip and is almost 8' from the aim-point, this 
long observation provides nearly 300 counts and allows the best 
study of the source spectrum. In this Letter, we concentrate on 
the analysis of these two latterly mentioned observations; we 
also make use of other short (~ 5ks) ACIS observations (see 
e.g., Di Stefano et al. 2002; Kong et al. 2002) to construct the 
long-term lightcurve. 

All ACIS data were telemetered in Faint mode, and were col- 
lected with frame transfer time of 3.2 s. In order to reduce the 
instrumental background, we screened the data to allow only 
photon energies in the range of 0.3-7 keV, and ASCA grades 
of 0, 2, 3, 4, and 6. In addition, ^ 300 s of data was ex- 
cluded from the 2001 October 5 observations due to increased 
background (likely due to increased particle flux from the so- 
lar wind). Data were reduced and analyzed with the Chandra 
X-ray Center CIAO v2.2 package ' and spectral analysis was 
performed by making use of XSPEC vl 1 .2 ^. Unless otherwise 
specified, all quoted errors are at the Icr level. 

3. ANALYSES AND RESULTS 



http://asc.harvard.edu/ciao/ 
^http://heas arc. gsfc.nasa.gov/docs/xanadu/xspec/index.html 
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3.1. X-ray Image 

CXOM31J004327.7+411829 was cleai'ly detected in both 
observations, but we determined the source position and mor- 
phology from the shorter 5ks image because the source was 
significantly closer to the aimpoint in this dataset. We exam- 
ined the aspect of this observation with the thread provided by 
the CXC ^; and find an offset of 1.34" and -1.43" in R.A. and 
Dec, respectively. After correcting for this aspect offset, the 
coordinates of the center of CXOM31 J004327.7H-41 1829 ai-e 
found to be R.A.=00'' 43"' 27.8', DEC.=H-41°18' 29" (J2000) 
with a positional error of about 0.5". The optical counterpart of 
CXOM31 J004327.7H-41 1829 is source 2-033 in Magnier et al. 
(1995), source 15 in d'Odorico et al. (1980) and source BA23 
in Blair et al. (1981). This object was identified as an irregu- 
lar, faint and high [SII]/Hq (~ 1) SNR. The optical extension 
is about 10"- 18", with a slight elongation along the DEC. di- 
rection. This source was previously detected (but not resolved) 
with the ROSAT PSPC and identified as an SNR (Supper et al. 
1997,2001, source RX J0043.4H-41 18). It was out of the field- 
of-view of the ROSAT HRl search for SNR in M31 which might 
have resolved it (Magnier et al. 1997). Figure 1 shows the "true 
color" X-ray image of CX0M31 J004327.7H-41 1829. The color 
scheme is defined as 0.5-0.8 keV = red, 0.8-1.2 keV = green, 
and 1.2-4 keV = blue. The choice of such energy ranges is to 
highlight the O and Ne/Fe emissions (see § 3.2) from the SNR. 
CXOM31 J004327.7H-411829 is clearly extended in X-rays as 
a distinct ring-shaped object with a diameter of ^ 11" (42 pc 
at 780 kpc). The uncertainty of physical dimension of the SNR 
is dominated by the distance to M31; for instance. Supper et 
al. (1997,2001) used 690 kpc, corresponding to ^ 37 pc. Our 
quoted distance (780 kpc; Stanek & Garnavich 1998) is the up- 
per end of all distance estimates. 
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Figure 1: "True color" Chandra AClS-1 image of 

CXOM31J004327.7-I-411829. This image was constructed from 
the soft (red: 0.5-0.8 keV), medium (green: 0.8-1.2 keV) and hard 
(blue: 1.2^ keV) energy bands. The pixel size is 0.496" and the 
image has been slightly smoothed with a 0.8"a Gaussian function. 
The diameter of the SNR is about 11". 



3.2. Spectral Analysis 

The long 37.7 ks ACIS observation was used for spectral 
analysis. We extracted data from a circle of 12 pixels 6") ra- 
dius centered on CXOM31 J004327.7H-41 1829 and background 
from an annulus with inner and outer radii of 30 and 60 pixels, 
respectively. We found 266 net source counts. In order to allow 
statistics to be used, the spectrum was grouped into at least 
15 counts per spectral bin. Response files were selected ac- 
cording to the CCD temperature with standard CI AO routines. 
We fit the data with several single-component spectral mod- 
els including power-law, thermal bremsstrahlung, blackbody, 
Raymond-Smith (RS) and non-equilibrium ionization (NEl) 
models with interstellar absorption. The RS model is a sim- 
ple coUisional equilibrium ionization model, while NEl model 
is appropriate for modeling SNRs whose age is smaller than 
the time required to reach ionization equilibrium. NEl model 
consists of electron temperature {kTe) and ionization timescale 
(r = riet), where vie and t are the mean electron density and the 
elapsed time after the plasma was shock heated to a constant 
temperature kTe- Both models are often applied to study X- 
ray emission from extragalactic SNRs (see e.g., Wang 1999; 
Schlegel, Blair, & Fesen 2000; Hughes, Hayashi, & Katsuji 
1998). We were also concerned that the quantum efficiency 
degradation of ACIS will affect our spectral fit"*. The degra- 
dation is shown to be a function of time and is most severe at 
low energies. We therefore applied corrections using the ACIS- 
ABS absorption model in XSPEC ^. This model allows us to 
correct the response by inputing the number of days between 
Chandra''^, launch and the observation (805 days for our obser- 
vation). In general, the degradation mainly affects the best fit 
Nh, which is over-estimated without the correction. All results 
reported here are corrected for the degradation. 

Except for the RS and NEl models, other spectral models 
give unacceptable fits to the data {y^ /v > 2). For the RS and 
NEl models, we first fixed the abundances at solar value (An- 
dens & Grevesse 1989). We also did fits with the abundances 
of O, N, and Fe as free parameters and the remaining elements 
fixed to the solar value. Finally, we fixed the abundances based 
on optical spectroscopy of CXOM31J004327.7H-41 1829 with 
to be 2.3 X 10-4 (O), 8.4 x 10"^ (N) and 0.72 x 10"^ (S) relative 
to hydrogen, respectively (Blair et al. 1982); these correspond 
to 0.27, 0.75, and 0.44 solar abundance. Table 1 lists the best- 
fitting parameters for the RS and NEl models. 

All models except the RS model fixed at the optically de- 
termined abundances provide acceptable fits. The NEl models 
clearly have relatively large errors, especially the electron tem- 
perature, kTe. The best-fit A^h is slightly lower (within la) than 
the Galactic value along the direction (7 x 10^** cm"-; Dickey 
& Lockman 1990), and the lower limit is not well constrained. 
Moreover, the abundance of O, Ne, and Fe is not constrained if 
they are left as free parameters in the NEl model. The best fit 
temperatures of the RS models with solar abundances and free 
O, Ne and Fe abundances are similar The Nh for the models 
with solar abundances is an order of magnitude greater than the 
optically determined value (see below). The X-ray spectrum 
with an absorbed RS model is shown in Figure 2. It is evident 
that the spectrum is dominated by a broad emission line of O 
VIll at 0.654 ke V and a blend of Fe L shell lines and Ne K shell 
lines around 0.9 keV. 



^htt p://asc. harvard.edu/mta/ASPECT/fix_offset/fix_offset.cgi 
'^see ittp://asc. harvard.edu/cal/Links/Acis/acis/Cal_prods/qeDeg/ 
^http;//www.as tro.psu.edu/users/chartas/xcontdir/xcont. html 
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Table 1 

Best-fitting Spectral Parameters 



Model Nh kTus/kT, 


logriet 






Abundance^ 




Lx " xi/d.o.f 


(xl02icm-2) (keV) 




N 


O 


Ne S 


Fe 





RS 6.78t^^ 0.14!0:03 

RS 4.60+0?^ 0.18+«[J| 

RS 5.88!«J0 O.U^f, 

NEI 0.2l!^-6f 2.1!5.6 10.3!jj;^ 

NEI 0.25t§ 2.1!5.6^ 10.3!«:j^ 



1 


1 


1 


1 


1 


20.7 


1.06/13 


1 




0.78!«i4 


1 




4.5 


0.89/10 


0.75'= 


0.27'= 


1 


0.44'= 


1 


13.4 


1.64/13 


1 


1 


1 


1 


1 


0.32 


0.83/12 


0.75'= 


0.27'= 


1 


0.44'= 


1 


0.32 


0.83/12 



Note. — All quoted uncertainties are la. ACISABS model is applied to correct the degradation of ACIS (see text). 
^Relative to the solar abundance. 

''0.3-7 keV luminosity (x 10^' erg s~^), assuming a distance of 780 kpc; Stanek & Gamavich 1998. 
'=Fixed at optical value (Blair et al. 1982). 
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Figure 2: Upper panel: The Chandra ACIS-S4 spectrum of 
CX0M31 J004327.7+411829 with an absorbed RS model {Nh = 4.6 x 

10^' cm"^, kTRs = 0.18 keV). Lower panel: Residuals after subtracting 
the fit from the data in units of counts s"' keV~' . 

The hydrogen absorption coluimi density for the renonant 
from optical observations was estimated to be A^^ = 7x10^" 
cm-2 (Av = 0.4; Blair et al. 1981, 1982). The best fit value for 
Nh from the RS (NEI) model is higher (slightly lower) than the 
optically determined value. We also fixed Nh at 7 x 10^" cm"^ 
and the fit is still acceptable, but gives a shghtiy higher thermal 
temperature for both models. 



3.3. Lightcurve 

The long-term lightcurve is shown in Figure 3, which is con- 
structed from a series of ACIS pointings (see e.g., Di Stefano et 
al. 2002; Kong et al. 2001). The background subtracted count 
rate was corrected for exposure, background variation and in- 
strumental PSF. The mean count rate is 0.007 count s"^ and it is 
evident that the source does not vary significantly as expected 
for a SNR. 



We also tested for variability over shorter time scales during 
the 37.7 ks observation. We binned the data into 500 s, 1000 s, 
5000 s and 10000 s resolution and used Kolmogorov-Smimov 
test to search for variability, finding nothing significant. 



0.01 



<u 0.008 
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Figure 3: Lightcurve of the Chandra count rate in the 0.3-7 keV 
range of CXOM31J004327.7+411829 since 2000 March. The hori- 
zontal line is the average count rate of the SNR. 

4. DISCUSSION 

CXOM31 J004327.7H-41 1829 is clearly resolved with Chan- 
dra into a ring-shaped object. This is the first extragalac- 
tic SNR resolved at X-ray wavelengths (other than those in 
the Large and Small Magellanic Clouds). The diameter of 
CXOM31 J004327.7+411829as measuredby Chandrais about 
11", corresponding to 42 pc; this is comparable to the optically 
determined value (d'Odoricoetal. 1980; Blair etal. 1981). The 
X-ray image shows evidence of a spatial dependence of chem- 
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ical composition, with O concentrating in the north-western 
part, and Ne/Fe distributing along the eastern and south-eastern 
arc. It is unclear if the X-ray image shows the "blow out" seen 
in the south-eastern section of the optical image (d'Odorico et 
al. 1980; Magineret al. 1995). There is no significant variabil- 
ity on both long (months to years) and short (hours) time scales. 
The X-ray emission from CXOM31J004327.7H-411829 can be 
acceptably fit with a RS model. For abundances fixed at solar 
values, the estimated luminosity is near the Eddington luminos- 
ity for a 1 AMq neutron star, which makes this SNR one of the 
most luminous (> lO-'^ erg s"') objects in M31 (see Kong et al. 
2002). Beacaue it is not realistic to assume solar abundances for 
all elements, we here consider the RS model with varing chem- 
ical abundances. Fits to RS model with variable abundances 
find a slightly lower X-ray luminosity of ^ 5 x lO"*^ ergs"' 
(0.3-7 keV; Nh = 4.6 x 10^' cm'^ and temperature kTus = 0.18 
keV), which is at the high end of luminosity distribution of SNR 
in M31 (Supper et al. 1997; Magnier et al. 1997). Assum- 
ing CXOM31 J004327.7H-411829 is in the adiabatic expansion 
phase, we can estimate its physical parameters through the Se- 
dov solution: 

R = 5pcx(E5itl/no)''\ (1) 

where R,Esi,t3, and no are the radius (in pc), initial explo- 
sion energy (in units of 10^' erg), age (in units of 1000 yr), and 
number density of ambient gas (in units of cm"-^) of the SNR, 
respectively. The shock temperature can be written as 

r, = (o.i8 keV)(/^/f3)^ (2) 

Adopting a radius of 21 pc, Ts = 0.18 keV, and £51 = 0.465 
(Blair et al. 1981); from equations (1) and (2), We obtain 
t = 21000 years and no = 0. 16 cm"^. Considering the SNR's size 
and uncertainty of spectral fit, f ranges between 15700 years 
and 22300 years while no varies between 0.11 cm"-* and 0.26 
cm-^ We therefore can classify CXOM31 J004327.7H-41 1829 
as a middle-aged SNR (cf e.g., Hughes et al. 1998) based 
on the age derived from the X-ray spectrum. Because SNR 
sometimes contain a pulsar and/or a pulsar nebula we also 
did some fits including a power-law component with photon 
index fixed at 2 along with our RS model. These fits limit 
the power-law component to less than about 25% of the flux. 
From the derived density, the swept-up mass is estimated to 
be (4/3)M//no7r/?-' ^ 153Mq, where Mh is mass of a hydro- 
gen atom. Taking into account of the uncertainty of the SNR's 
size and temperature, the swept-up mass ranges from 105Mq 
to 172Mq. Thus for an ejected mass of a few solar masses, the 
ratio of swept-up to ejected mass is ^ 20-30, indicating that 
the SNR is in the Sedov phase. 



The X-ray emission from CXOM31 J004327.7H-41 1829 can 
also be fit with an absorbed NEI model although the uncertain- 
ties of Nh and kT^ are large. Using this model, the X-ray lu- 
minosity (0.3-7 keV) of the SNR is estimated to be ~ 3 x 10^'' 
ergs"', which is typical of SNR in M31. Using Eqns (1) and 
(2), we obtain t = 3210-7530 yeai's and no = 0.003-0.02 cm'^ 
The ratio of swept-up to ejected mass is ^ 5, indicating that 
the SNR might be in the early Sedov phase. Comparing to RS 
model, these values are smaller. In particular, the number den- 
sity of ambient gas is very low. The relatively high kT^ (=2.1 
keV) is not unusual; for example, N132D in the Large Mag- 
ellanic Cloud requires a two-component NEI model for which 
the temperatures are 0.8 and 2.7 keV (Favata et al. 1997). The 
usual explanation for the hot component is that it comes from 
the shock-heated swept-up circumstellar medium, or it is due to 
the inhomogeneity of the ISM. 

The number density of ambient gas is relatively low (no = 
0.003-0.26 cm"^). For example, from a well-studied sample 
of SNRs in the Large Magellanic Cloud (e.g. Hughes et al. 
1998), no ranges between 0.3 to 4 cm"^. Previous ROSAT oh- 
servations of M3 1 indicate that the number of X-ray emitting 
SNRs is surprisingly small (Supper et al. 1997; Magnier et al. 
1997); only 16 out of 200 opticaUy identified SNRs were de- 
tected in X-rays. Magnier et al. (1997) suggested that this lack 
of X-ray emission could be due to lower than expected local 
ISM densities, i.e., non-detected SNRs may be in regions with 
densities less than 0.1 cm"^^, while those SNR detected in X- 
rays may be in higher density regions. The low no we find for 
CXOM31J004327.7H-411829 provides support for this sugges- 
tion. In addition, in the optical image of Magnier et al. (1995), 
the SNR appears to be away from high density regions of the 
galaxy, suggesting that it may be in a relatively low density re- 
gion. 

Magnier et al. (1997) note that a multi-component ISM 
model (McKee & Ostriker 1977) is required for M31. In this 
case the densities we derive herein may be upper limits to the 
dominant low-density component of the ISM in the vicinity of 
the SNR. A future longer observation with the SNR on-axis 
is required to investigate multi-component model and detailed 
morphology. 
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